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Abstract 

Quantum mechanics abounds in models with Hamiltonian operators which are 
energy-dependent. A linearization of the underlying Schrodinger equation with 
H = H{E) is proposed here via an introduction of a doublet of separate energy- 
independent representatives K and L of the respective right and left action of 
H{E). Both these new operators are non-Hermitian so that our formalism ad- 
mits a natural extension to non-Hermitian initial H{E)s. Its applicability may 
range from pragmatic phenomenology and variational calculations (where all 
the subspace-projected effective operators depend on energy by construction) 
up to perturbation theory and quasi-exact constructions. 
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1 Introduction 



In non-relativistic quantum mechanics the role of time is exceptional as it 
remains, up to a few really exotic models [1], unquantized. This makes the 
current Schrodinger equation i % dt'4>{t) — H ■4>{t) easily integrated in all the 
time-independent cases where we get ipit) = exp[i {t — to) H] ipito). This means 
that we just have to solve the time-independent Schrodinger equation 

H I = I (1) 

and evaluate using the spectral representation of the Hamiltonian opera- 
tor, 

= ^ (2) 

n 

The picture becomes less transparent in all the time- dependent regimes where 
we have to deal, e.g., with the time- dependent boundary conditions [2] and 
where the definition of H = H{t) becomes transformed and transferred to 
H — H[E) in an appropriate Fourier-transformation modification of eq. (1), 

H (e^")) |^(")) = |^(")) . (3) 

In our present short note we intend to make a few comments concerning the 
latter "non-linear" Schrodinger equation. 

2 Motivation 

2.1 An emergence of the energy-dependence in physics 

The energy-dependent version (3) of the current bound-state Schrodinger equa- 
tion (1) is far from being an exotic problem without any practical relevance. 
Almost an opposite statement is true, at least in principle. Indeed, whenever we 
contemplate a "real" system and construct its manifestly energy-independent 
"realistic" Hamiltonian Hr with due care, we must keep trace of many "possi- 
bly relevant" degrees of freedom. In the next step we usually decide, on some 
more or less intuitive grounds, that the "really relevant" part of might be 
just its projection H(^ms) = P Hr P on a certain "model" Hilbert subspace. 
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Having finished such a procedure (which is quite common, say, in nuclear 
physics [3]) one is obhged to compare the spectra {Eji} (of Hji) and {E(^ms)} 
(of Hi^MS)) and to verify, in this way, the quahty and/or rehabihty of our initial 
intuition. A formal tool is provided by perturbation theory [4] in terms of 
which we may arrive at a more quantitative verification of our hypotheses by 
the almost elementary replacement of the "exact" equation Hr \^ r) — Er \^ r) 
by its two projections P and Q — 1 — P, 

P {Hr - Er) P I^k) + P{Hr- Er) Q \^r) = 0, 
Q (Hr-Er) P\^r) + Q (Hr-Er) Q\^r)^0. 

Eliminating the irrelevant projection Ql^/?) we easily arrive at the strictly 
equivalent reduced eigenvalue problem H(^red){ER) \ (p) — Er \ if) of the form (3) 
with 

Q 



H[red){E) = < Hr + HrQ 



QHr). (4) 



Q {Hr - E) Q 

Using such a construction of H(^red) = H(^red){E), all the energy errors vanish 
and we have just the identical energy spectra {Er} = {E(^red)}, in spite of the 
arbitrary reduction \^r) — > P \^r) = \(f) oi the Hilbert space itself. 



2.2 Energy- dependent models in computations 

The role of the above exact formula (4) in practical computations is indirect 
because the energy-dependent operator difference between the exact Hr and 
reduced H(^red)iER) is, undoubtedly, hardly obtainable in closed form. This 
means that the unique and exact parameter-dependence of H(^red){^) is too com- 
plicated. In practical computations it invariably requires a perturbative, vari- 
ational or at least empirical or trial-and-error simplification H(^red){^) H{z). 
At all the relevant real values of z this re-opens the possibility of employing 
the spectral decomposition, 

H{z) = E 1^^"^^)) E^-\z) {¥-\z)\ , G {4™^}~^„ • (5) 

n 

We see that the work with the energy-dependent Hamiltonians may be both 
feasible and useful [5]. One only has to keep in mind that it still requires 
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comparatively lengthy calculations because at any pre-selected superscript we 
must guarantee that the value of z satisfies the self-consistency condition 

^phys — -S'^ ^-^-p/ij/s) • (6) 

This means that the algorithm of construction of bound states generated by 
any energy-dependent Hamiltonian H{E) remains non-linear, with all the un- 
pleasant mathematical consequences including the possible non-existence or 
redundancy of the real (and even complex) solutions of eq. (6) at every n. 

3 Examples of applicability 

3.1 Energy scale and its subdomains 
3.1.1 Thresholds in phenomenological physics 

In the majority of phenomenological models Hp^,, one does not specify any 
upper limit E^ax of their applicability. Thus, atomic and nuclear systems are 
very often described by non-relativistic Schrodinger equations with interactions 

reduced to a local two-body phenomenological potential V = V{xi — Xj) [6]. Of 
course, the credibility of such a model becomes very low at the higher energies 
so that we tacitly have to assume the existence of a complete, two-component 
description 

where the precise localization of the onset of the new physics (expected to occur 
somewhere near the "threshold" energy Ecru. ~ 300Me\/ in nuclear physics) 
is rarely considered really important because, for years, the detailed structure 
of Hhighen. has been believed to be prohibitively complicated. Still, due to the 
steady increase of our computing power, the high-energy parts of eq. (7) and 
of its step-function energy-dependence generalizations 

H = Ha,bioTEe{Ea,Eh) (8) 

might acquire the status of a tractable problem quite soon [5, 7]. 
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3.1.2 Smoothly dependent Hamiltonians 

One of the best known theoretical sources of differences between the respec- 
tive low- and high-energy Hamiltonians are relativistic phenomena (e.g., the 
decays and the emergence of various new degrees of freedom) and corrections 
(say, in the form of a weakly energy-dependent effective mass in nuclear [8] and 
particle [9] physics). In such a setting, the building of a model should start 
in the relativistic kinematical regime and proceed towards the lower energies 
by a systematic simplification of Hhighen.- One usually works with a well de- 
fined energy-dependence in Hfiighen {E) which is smooth. A wealth of explicit 
examples of this type fits in the scheme of eq. (8) and occurs in numerous 
applications [10]. 

Naturally, the study of the models H{E) with a smooth energy dependence 
attracts attention also by its methodic aspects. To grasp their flavor, the reader 
is recommended to play with the harmonic-oscillator example of ref . [7] , 

(9) 

an extremely easy solution of which illustrates many paradoxes exhibited by 
much more complicated realistic models. 

3.2 Energy-dependence as a formal freedom 
3.2.1 Non-linectrity and non-orthogonality 

The key constraint (6) may be satisfied at a large as well as empty set of real 
roots at a given H{E) and n, 

^^^^^ _ ^(n,l)^ ^(n,2)^ ^ ^(n,m(n)) _ 

In a compactified notation using multi-indices a — {n,j) we shall abbreviate 
^{n,j) _ denote \^^'^\zphys)) — |0a) in order to identify characteristic 

difficulties which arise in connection with the energy dependence in if = H{E). 

First of all, we have to emphasize that the standard orthogonality relations 
between the separate bound states become lost. Even though the norm ||0a|| 
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of each particular eigenstate may easily be fixed through evaluation of the self- 
overlaps (0a 1 0a), it is necessary to evaluate also all the off-diagonal overlaps 

(0a|0/3) = (10) 

which need not non- vanish in general. Secondly, whenever we have to deal with 
an energy-dependent model where the energies remain discrete, we may alter 
the denotation ol H = H{E^) = H{E^'''^^) = i/(")(^("'-')). In particular, the 
latter convention may prove suitable whenever the partial energy-independence 
property 

H{n)(^^{n,j)^^ jj{n)(^^{n,k)^^ J, = 1, 2, . . . , m(n) (11) 

is encountered as a generalization of the step-shaped energy dependence (8). 
3.2.2 Quasi-exact solvability 

We shall see below that one of the most important technical assumptions of the 
feasibility of the work with any H = H{E) is the feasibility of the evaluation 
of the overlap matrix R (10). This is only trivial in the completely solvable 
and fully energy-independent limit of eq. (9) (with m{E) = const and with the 
safely diagonal -Ra,/?) where the evaluation of the sequence of the non- vanishing 
overlaps may be reduced to the formula 

j^iharm.oscil.) ^ ^ c„ / C"^' X^'^^^+^^dx (12) 
n=0 "''^ 

which evaluates to a sum of F— functions in both the Hermitian [11] and non- 
Hermitian [12] cases. Difficulties with the determination of R perceivably 
increase whenever one tries to move to a more realistic model. Vice versa, 
the requirement of the preservation of a manifestly non-numerical form of the 
overlaps R leads directly to one of the eligible definitions of the concept of 
the so called quasi-exact (QE) solvability [13]. For an exphcit illustration of 
this statement, let us pick up one of the most popular QE models, viz., the 
sextic QE anharmonic-oscillator generalization of eq. (9) which, in the units 
h'^ = 2m{E) = C = 1, reads 
d2 



--^+A(Ea) r^ + Br^' + Cr'' 
dr^ 



*a(r) = £;a*a(r) • (13) 
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This QE model was discovered by Singh et al as early as in the late seventies 
[14] but its analytic continuation to a non-Hermitian regime appeared only 
very recently [15]. In both these cases we may use our above composite-index 
notation with a = {N,j) and j = 1,2, ... ,m{N) and emphasize that the 
bound-state solutions of eq. (13) remain "quasi-exact" (i.e., elementary and 
proportional to polynomials) if and only if we keep the integer fixed. We 
must select an appropriate specific "spring constant" A{Ea) — A{E^^'^^) = An 
which is different for the different N — 1,2, . . .. In this manner one introduces 
a certain "minimally complicated" energy dependence in the Hamiltonian (13). 
At the same time, all the sums and integrals which define the diagonal as well 
as non-diagonal overlaps i?^,/? remain very similar to eq. (12) and may be easily 
shown to degenerate to finite sums of incomplete F— functions. 

4 Hermitian Hamiltonians H(E) 
4.1 Bi-orthogonality and completeness 

In eq. (10), a nontrivial assumption is needed as a guarantee that the matrix 
R has an inverse (for the time being, let us also keep in mind that it is, by 
construction and/or assumption, Hermitian). Only then, the formal decom- 
position of the identity projector becomes available as a double sum over the 
suitable (sub)set A of the multi-indices. 



Even if we knew the set A and proved the convergence, the practical value of 
such a formula would remain significantly lowered by the non-diagonality of 
the infinite-dimensional matrix R. The sufficiently precise evaluation of this 

matrix and of its inverse R~^ is needed. As we already mentioned (cf. also 
[16]), all this requires a sufficiently elementary form of the basis states 

Equation (14) acquires the standard linear- algebraic meaning of a "com- 
pleteness" relation for our (non-orthogonal) basis { }aeA- After we abbre- 




(14) 
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viate 

l3eA 

a new basis emerges which, by construction, obeys the Kronecker-delta-overlap 
rule 

{{(t>a\^l3) = {(t>a\(t>l3)) = (^a,/?, a, (3 & A 

and is called, on this ground, bi-orthogonal and bi-orthonormal [17]. 
4.2 Two alternative quasi- Hermitian representants 

Practical use of the bi-orthogonally generalized basis exhibits a lot of parallels 
with its ordinary orthogonal predecessor. First of all, it enables us to re- 
interpret the double series (14) as a single-index expansion employing the two 
different types of vectors, 

\4>a){{4>a\ = Y. \'Pa)){K\- 

It is easy to verify that the new operator defined by its generalized spectral 
representation 

or, in an abbreviated notation, 

K^Y. \4>a)E^{{4>^\ 
aeA 

may be re-interpreted as an operator which shares with H{E) its action to the 
right, 

K\(t)^)^E^\(t>^) , aeA. 

By construction, the new operator K is energy-independent (i.e., no nonlincar- 
ity is encountered). The latter property is counterbalanced by non-Hermiticity 
K ^ which means that the action of K to the left is much more compli- 
cated. Nevertheless, we may complement the linear, non-Hermitian operator 
K by another operator, viz., by 

a,/3eA ''^ 
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This means that 

l3eA 

so that the second non-Hermitian hnear operator L shares with H{E) its action 
to the left. In this sense, the single nonhnear operator H{E) may be understood 

as represented by the linear doublet {K,L). 

One has to notice that K = and L = K'^ are mutually inter-related, 

or, alternatively, 

Lr'=r'i^= E I0a»^a((<^a|, = E l'^«)>((0a|- (16) 

As long as wc have ,^ > 0, the terminology of the review paper [18] may be 
used implying that both our mutually conjugate Hamilton-like operators K 
and L are quasi-Hermitian. For this reason, each of them admits a consistent 
quantum-mechanical interpretation (readers are recommended to check ref. 
[18] for more details). 



5 Non-Hermitian Hamiltonians H{E) 

Recently, several groups of authors [19] - [24] tried to weaken the traditional 
Hermiticity H = [tacitly also assumed, up to now, in formulae (2) and (5) 
above] . The standard pattern of such a generalization (which, for the sake of 
simplicity, keeps the energies real) lies in the assumption that the bra and ket 
vectors in the similar spectral expansions are not the mere Hermitian conju- 
gates of each other. This resembles the above-described relation between the 
vectors [(^q) and |0a)) which we also assumed to remain significantly different. 
Indeed, many rules of our preceding section 4 will find their parallels also in 
the forthcoming text where we shall admit that H[E) ^ W{E). 

5.1 Quasi-Hermitian input Hamiltonians H{E) 

Let us keep the supers criptedkets l^^")) of eq. (2) unchanged and introduce an- 
other, independent infinite series of their subscripted paLitnevs {"^(n)) 7^ |\1''-"-'). 
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In the light of the current textbooks [17] the usual orthogonality (and normal- 
ization) assumption (^("*)|^(")) = 5jn,n finds its most natural generalization in 
the so called bi-orthogonality (or rather bi-orthonormality) relations between 

these two sets, 

= m,n= 1,2,... . 

In parallel, the usual completeness relations for a basis in our Hilbert space 
Ti. must be replaced by the following innovated, bi-orthogonal-basis-related 
formula 

/ = ^|v[/W)(x[/(„)|. (17) 

n 

This enables us to replace the spectral decomposition (2) of Hermitian H — 
by its non-Hermitian bi-orthogonal-basis analogue or generalization, 

^^ = E (18) 

n 

Parallels are preserved with eqs. (16) and/or (15), 

n n 

Hrj-^ = r/-^ = E 1*^"^) E^""^ = E l*^''^) • (20) 

n n 

We may call again all our generalized Hamiltonians (18) quasi- Hermitian [18] 
since we are working with the manifestly regular and positive > in the 
latter two relations. 

5.2 Innovated pair of the represent ants K and L 

Once wc keep in mind the quasi-Hermiticity symmetries (19) and (20), wc 
may introduce an error-checking convention under which we only consider the 
formulae where all the kets are upper-indexed while all the bras are lower- 
indexed. Under this convention the extension of the results of section 4 to all 
the quasi-Hermitian H{E) ^ H.\E') is straightforward. Firstly, in the way 
inspired by eq. (5) we have to define 

i/(z)=^]*W(z))E(")(z)(*(„)(z)| (21) 

n 
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and re-abbreviate 



W'^) = \^!^^\z,nys)), {Va\ = (*(n)(%fe,.)| (22) 

where the rising or lowering of the index means a transition to another, entirely 
different vector. Assuming the knowledge of all the overlaps 

(notice an inessential change in our matrix indexing), we have to emphasize that 
they do not form a Hermitian matrix anymore. Still, its assumed invertibility 
suffices for us to write down the unit projector 

i = E (^"0/ (^/^l = E W')) (v^/^l = E Iv'") ((^al (23) 

where we abbreviated 

aeA l3eA 

This means that 



and that 



^=E l^niM^i: 

aeA aeA 



It is now trivial to define the generalized operators K and L which share with 
the quasi- Hermitian H{E) its action to the right and left, respectively, 

aeA 

peA 

This completes our construction. 
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6 Summary: Separation of the left and right 
action of the energy-dependent Hamiltoni- 
ans in Hilbert space 

When we compare the results of sections 4 and 5 we may feel surprised by 
the smoothness of the transition to the non-Hermitian H{E). In fact, the 
only perceivable consequence of the loss of the Hermiticity of H{E) lies in the 
emergence of the independent operator of the inverted metric which is defined 
in terms of the new, upper-indexed bras and kets. Formally, the breakdown 
K ^ and L ^ of the conjugation symmetry only implies that the 
quasi- Hermiticity rules (19) and (20) must be replaced by their appropriate 
modifications 



and 



aeA aeA 



where we must keep in mind that 

- E 1^") ' = E 1^")) ((^1 • (24) 

aeA aeA 

This means that we in fact violated the "error-correcting" convention of section 
5 and started writing all the "redundant" Hermitian-conjugate formulae in 
their explicit form. Such a step facilitates the most concise formulation of 
our present message stating that the original operator H{E) and its conjugate 
H\E) are in fact subtly different in both the Hermitian and non-Hermitian 
cases. In the other words, even after the present "linearization" of their action 
in the standard, self-dual Hilbert space, none of their two energy-independent 
representatives K and L is redundant. 
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